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We revisit the quantum lattice gas model of a spinor quantum field theory — the smallest scale 
particle dynamics is partitioned into unitary collide and stream operations. The construction is 
covariant (on all scales down to a small length I and small time t = cl) with respect to Lorentz 
transformations. The mass m and momentum p of the modeled Dirac particle depend on I according 
to newfound relations m — ttio cos 2^ and p — jsin2|l, respectively, where A is the Compton 
wavelength of the modeled particle. These relations represent departures from a relativistically 
invariant mass and the de Broglie relation — when taken as quantifying numerical errors the model 
is physically accurate when ^ ^ A. Calculating the vacuum energy in the special case of a massless 
spinor field, we find that it vanishes (or can have a small positive value) for a sufficiently large wave 
number cutoff. This is a marked departure from the usual behavior of such a massless field. 
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I. INTRODUCTION 



We consider the dynamics of a spinor quantum field 
where spacetime becomes discrete at scales smaller than 
some fundamental length. In particular, we revisit the 
quantum computational lattice representation known as 
the quantum lattice gas model, a dynamical Feynman 
chessboard model of the Dirac equation Variations, 
rediscoveries and improvements of the Feynman chess- 
board model have appeared over the years, including a 
model in 3+1 dimensions by Riazanov an Ising spin 
chain representation by Jacobson and Schulmanlj, |5| , a 
fundamental deterministic model by 't Hooft 017 a lat- 
tice Boltzmann model by Succi and Benzi quantum Q, 
a unitary model by Bialynicki-Birula and quan tum 
lattice gas models in 1+1 dimensions by Meyer [l3| and 
in 3+1 dimensions by this author [ll[. We presently con- 
sider a representation that retains 4-momentum conser- 
vation E"^ = {cpf' + {mc^Y' of special relativity down to 
a small length scale t and time scale r. The low-energy 
limit of the lattice model is defined as the dynamical 
regime where the Compton wavelength A of the quantum 
particle represented by an amplitude field ipi^) is much 
larger than the small scale, ^{x) is treated as continuous 
for A ^ £. Continuous derivatives emerge as effective 
quantum operators and the particle physics may be repre- 
sented by the Dirac Lagrangian /^oirac = i^{i^^d^i — 'mo)4', 
where the Dirac matrices are 7'' = (7", 7'), the spacetime 
derivatives are d^^ = (9t, di) for i = 1, 2, 3, and rUo is the 
"invariant" particle mass (here expressed in natural units 
with h = 1 and c = 1 for convenience). 

This paper is organized as follows. We begin in Sec. [TT] 
by formally introducing the quantum lattice gas model as 
a Langanian based theory. Then, in Sec. IIIIl we present 
a mapping procedure whereby the discrete dynamics of 
a quantum lattice gas model is made equivalent to the 
Dirac equation. This procedure leads to analytical form 
of the particle momentum that is a modification of the 



de Broglie relation. In Sec. IIVI we present a deviation of 
the quantum lattice gas stream and collide operators that 
form the basis of our quantum algorithm for the Dirac 
equation. In particular, we derive a unitary expression 
for the collision operator that is serves as a mechanism 
to give the spinor field its mass. Our primary intent is to 
show that the quantum lattice gas, taken as a numerical 
tool for this quantum computational physics application, 
provides a high degree of numerical accuracy. Then, in 
Sec. [Vj we examine the newfound requirements to have 
the dynamical equation of motion of the quantum lattice 
gas model equal the Dirac equation at a selected small 
scale and explore the consequences of these requirements. 
We present a calculation of the vacuum energy density 
of a spinor field, treating the special case of a massless 
spinor field. Following a detailed analysis of the behavior 
the error terms, one finds an alternate theoretical purpose 
of the quantum lattice g toy model. It provides an 

example where the vacuum energy of a massless spinor 
field can vanish or be very small. That is, one can take 
the quantum lattice g toy model of Planckian scale 

physics and thus set t he sma ll scale sizes ^ and r to the 
Planck length = ^JhGjc? and Planck time Tp = ^/c, 
providing a route for a small positive cosmological con- 
stant. In Sec. IVIl is a conclusion and summary. 



II. QUANTUM LATTICE GAS MODEL 

The proposed high-energy (small scale) quantum lat- 
tice gas representation may be formally expressed by the 
Lagrangian density of the form 
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V' = /:o™.+o(t2). 

(la) 
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By the least action principle, this Lagrangian density 
leads to the equation of motion of the form 



-T7 



(lb) 



Equation (jlbp is the equation of motion of a quantum 
lattice gas, a unitary model for a system of noninter- 
acting Dirac particles. On the right-hand side of (fTbl . 
free chiral particle motion is given by a stream operator 



, with momentum operator pi 



-id,. A 



mass-generating interaction that causes a lefthanded par- 
ticle to flip into a right-handed particle (and vice versa) is 
given by a unitary collision operator Uc = e"''^™'*' . The 
product C/qlg = UsUc is the local evolution operator of 
a quantum lattice gas system acting on the spinor field 
^T(^) = (-^^(3.) iIj^(x)). The lefthand side of (HbJ is a 
newly computed state ip'{x) = e"^^* ^(2;) at time t + r, so 



lb]) may be written as a quantum algorithmic map 



(Ic) 



taken to be homogeneously applied at all points x of 
space and at all increments t of time. In natural units 
{h = 1 and c = 1), the quantum lattice gas model ([Tc]) 
is specified in 1 -I- 1 dimensions by the following unitary 
operators: 



Ua 



(2a) 
(2b) 



where mo is the mass of the modeled Dirac particle in 
the low-energy limit £/X ~ 0. 

In the low-energy limit, the C(t^) error terms on the 
righthand side of pa]) become negligible, so L'^"'^ ^ -Coirac 
is covariant with respect to Lorentz transformations in 
this limit. Yet, can £'^"'' be manifestly covariant at high- 
energies £/X ~ 1? We consider how to achieve the co- 
variance of ([Ta)) at a small scale: it necessarily occurs 
when the high-energy equation (jlbp — or equivalently the 
quantum lattice gas equation Px)) — has the form of the 
Dirac equation {"^^Pfj, + m)'ip{x) = 0. 

The model is an expression of the simple idea of a 
spacetime manifold that becomes discrete below a small 
scale i. The prescriptions needed to make px)) equivalent 
to the Dirac equation are derived in the following section. 



rotations" Ua^ 



and Uf, 



-^-^^ (3a) 

cos — cos — - — sm — sm —n\ ■ n-i 
2 2 2 2 

sm — cos —1X1+ cos — sm -^n^ 



. Pi . P2. 

sm — sm "2"'^! ^ 



(T, 



(3b) 



where cr = {a'x,<^y,<^z) is a vector of Pauli matrices, ni 
and 77.2 are unit vectors specifying the respective princi- 
pal axes of rotation, and /3i and /32 are real-valued rota- 
tion angles. Let us take = e ' 2 "-a cr quj- stream 

operator and Uc ~ e * 2 "i "^ as our collision operator. 
Without loss of generality, we may choose the principle 
axis of rotation along the z to generate streaming. 



■ ^2 



(4a) 



and treat the quantum algorithmic map as if it were ap- 
plied in 1-f 1 dimensions.^ In this frame a general collision 
operator is 



Vc 



i ^ (ao-i: +/3cr„ +70-2 ) 



(4b) 



where a, /3, and 7 are real valued components subject 
to the constraint + -f- 7^ ~ 1 . The unitary opera- 
tors (ID) are applied locally and homogeneously at all the 
points in the system. That is, we consider a construc- 
tion whereby the two principal unit vectors specifying 
the axes of rotation are 



fl2 = (0,0,1). 



(5) 



With this choice, nixn2 = (/3, —a, 0) and rii •n2 = 7, so 
([3]) is a quite general representation of a quantum lattice 
gas evolution operator 



U^Uc 



cos COS 

2 2 



. Pi . h 
7 sm — sm — 
' 2 2 



(6) 



. Pi P2 
a sm — COS — 
2 2 

„ . Pi P2 
p sm — cos — 
'^2 2 

. Pi P2 
7 sm — cos — 
' 2 2 



„ . Pi . P2 
' p sm — sm — 
'^2 2 

. Pi . P2 
' a sm — sm — 
2 2 

Pi . P2\ 

COS — sm — 
2 2 



III. IMPOSING COVARIANCE AT THE SMALL 
SCALE 

Here we show that (|lcp in the high-energy limit can 
be made equivalent to the Dirac equation. We begin 
with a local evolution operator as a composition of "qubit 



^ In l l3bl l wo used the identity (a • <t) • (b • cr) = ab + i (ax b) ■ cr. 

^ The reduction from 3-|-l to 1-1-1 dimensions is allowed because 
the algorithm has the product form tj}' {x) = UsUctpi^) 
where Us = e-'f '^I'C/f e'f (^!'+'^-'!7|e"'T'''-!7| = e"-^" ■"'"Pi , 
with Dirac matrices 7" = cr^. ig) 1 and 7* = icTy (g) cr^ in the chiral 
representation [3. 1 1 ill . Streaming in each of the spatial directions 
occurs independently, so for simplicity we can choose to consider 
a Dirac wave moving along z. 
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The Dirac equation for a relativistic quantum particle 
of mass Too may be written as 

ihdtip = -cpzazip + moC^cTxip. (7) 

Its time- difference form may be written as 



1 H Z CTz Z CTx ] V, (8) 



for small r and for momentum operator ~ —ihdz- We 
may view the unitary operator acting on the right-hand 
side of (ISl) as the effective low-energy operator obtained 
from the quantum lattice gas operator (jG)) 



small £ , ICPzT nJloC T 
UsUc >l + ——(7z 7 cr^. (9) 



To establish a correspondence between ^ and we 
simply choose the real- valued components of rii to satisfy 
the following three conditions: 



(10a) 
(10b) 
(10c) 



a sin 
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7 sin 
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+ COS 
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Additionally, we should respect the reality condition that 
fii have unit norm^ 



(lOd) 



that we established above with the collision operator 
(|4b[) . For the sake of simplicity, let us start with a spe- 
cialized construction whereby ni is perpendicular to n2 . 
The solution of (|10l) in this special case is 

a = cos— = — sm— 7 = 0. (11) 



Inserting ([TT|) into (|10ap gives sin ^ 



turn una) is ^/l - ( I sm^2 



and in 



In giD we 



chose —£pz/h = /?2/2, so in turn we have 



1 - 



moC^rN . ipz cpzT 



(12) 



in an intrinsic way. Equation (|12|) can be interpreted 
as a rather fundamental relativistic relationship between 
particles and points. In place of the theory of special 
relativity for classical particle dynamics in a continuum, 
here we have constructed a lattice-based version of special 
relativity for particle dynamics emerging at a small scale 
where the spacetime foam has a regular structure. 

Let us consider some implications of (HH. Squaring 
(IT2]) gives 



n . ipz 
— sm — — 
T h 



sin 



^Pz 

h 



^ ={cpzY. (13a) 



Then adding TOqC"* to both sides, we have 



■ 2 ^Vz ( 2 

sm — — h (TOqC 



,2\2 



{cpz) + (TOoC ) 



(13b) 

This is a candidate grid-level relativistic energy equation 
that leads us to define a grid momentum pl'"* and a grid 
mass TO dependent on i as follows: 



h . IPz 
— sm — — 
cr h 



TO = Too cos (14) 



^Pz_ 

h 



Hence, the lefthand side of (|13bp can be interpreted as a 
redefinition of the Dirac particle's kinetic and rest ener- 
gies. Inserting the de Broglie relation (p^ = hjX momen- 
tum eigenvalue), the grid mass and momentum become 



h . 2Tr£ 
— sm — — 
cr A 



27Te 



(15) 



In the low-energy limit defined by A ^ 2tt£, expanding 
([T^ to first order implies that the space and time cell 
sizes are linearly related by the speed light i = ct, an 
intuitive relationship that we expect to hold. In the low- 
energy limit, (|15p reduces to 



h 
A 



(16) 



That is, the low-energy limit of (|15p corresponds to a 
usual quantum particle with an invariant mass that is 
entirely independent of the particle's momentum, and 
the quantum particle acts like a wave according to stan- 
dard quantum mechanics. However, there is a marked 
departure from standard quantum mechanics in the high- 
energy limit in the region A < 20£ as shown in Fig. [1] 



This is a grid equation that relates the cell sizes £ and 
T to the mass and momentum of the quantum particle 



^ Alternatively, instead of ijlOdll . we could impose the condition 
that cos ^ cos ^ — 7 sin ^ sin ^ = 1, forcing ^ to be identical 
to lO. However, in this case, the resulting solution for compo- 
nents of ni has a imaginary, and this breaks the unitarity of Uc- 
So, we impose l llOdl l to strictly enforce unitarity. 



IV. THE ALGORITHM IN NATURAL UNITS 

For expediency, let us now switch our dimensional con- 
vention to the natural units, h = 1 and c = 1.^ We can 



* In the natural units h = 1 and c 
like dimension of length (i.e. [£] = [- 



1, length and time have 
= L) while mass, momen- 
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write (IT2I) as 



a/I — TOgT^ C0S( 



(17a) 

,) = v/l-^V2, (17b) 
Fur- 



or equivalently as e = exp 

thermore, our solution (fTTj) implies that the rotation axis 
for the collision operator is ri-i = x cos ipz + y sin (.pz , and 
in turn this implies that the hermitian generator of ()4bp 
is ni • cr = a ^ cos Ipz + <7ySm£pz = axe'^"'' . Hence, 

since sin ^ = TOoT, we can explicitly represent the colli- 

■ til - 

sion operator Uc ^ e * 2 "1'^ in terms of the mass and 
momentum of the quantum particle as 



Uc 



(18a) 
= moT. (18b) 



Multiplying by the stream operator C/| 
evolution operator ([6]) can now be explicitly calculated 



ie'^P^'^^cr^e^'^-^f^ mor(19a) 



(19b) 



■\/l — mfr^ cosp^f + i(T;j-\/l — m^r^ sinp^f 



JTTbl 



This result leads us to define the rotation axis 
ni2 



mo . Pz . 

X -\ z. 

E E 



(19c) 

(19d) 
(19e) 



(20) 



Since (ni2 • cr) = 1 (an involution), we are free to write 
5e]) in a manifestly unitary form e ' 2 "•12-'^ a,s follows: 



exp 



Vl ) ni2 • cr 



exp 



cos ^ \/l — E'^t'^ 
I [o-zPz - Oxmo) 



E 

where in the last line we made the identification 



cos(i;£) = Vl - E^T^, 



(21a) 
(21b) 
(21c) 

(22) 



turn, and energy values have like dimension of inverse length (i.e. 
[m] = [p] = [E] = L~^). Any expression written in the natural 
units can be converted back to an expression in the dimension- 
ful M,L,T units by simply reinserting the speed of light and 
Planck's constant by £ i->- | , m 1-^ 2^ ^ p ££ ^ and E 1-^ ^. 



which is exact to one part in 10^^^ (i.e. accurate to 4th 
order in £). This is equivalent to identifying the gate an- 
gle with the ratio of the small scale length to the particle 
energy, -^2/2 = i/E. Then, since U^Uc = e"*''""''^, 
the high-energy Hamiltonian may be written as 



- I -PzO-z + mocr.x 



(23) 



Thus, we have successfully demonstrated that the quan- 
tum lattice gas equation (jlcp is equivalent to the Dirac 
equation as the hamiltonian generating its unitary dy- 
namics is the Dirac hamiltonian, even at the small scale. 



V. CONSEQUENCES OF THE MODIFIED DE 
BROGLIE RELATION 

Imposing covariance on the high-energy representation 
([T|) leads to two departures ([15]), that we derived in the 
previous section, from the correct behavior given by rel- 
ativistic quantum field theory. First, the momentum of 
the quantum particle must obey a small length £ depen- 
dent momentum relation 



h . 2n£ 
p = — sin ■ 



£ 



A 



(24a) 



in place of the de Broglie relation p ~ h/\. Second, 
the mass of the Dirac particle is no longer taken as an 
invariant quantity — it must depend on the small length 
as well as the particle's Compton wavelength 



m = mo cos ■ 



27r£ 



(24b) 



where mo is a fixed constant, otherwise interpreted in the 
low-energy limit as the invariant particle mass. Plots of 
are given in Fig. [T] for mo set to the electron mass. 
Notice that (|24ap vanishes for A = ^ and A = 2^ and 
oscillates about zero as A — > 0. Also notice that (|24b[) 
vanishes at A = 4^ and is in fact negative for A = 2^ and 
A = 3£, again oscillating as A — >■ 0. These departures 
from standard quantum mechanics and special relativity 
theory are rather consequential at very small scales, de- 
parting on scales < 2Q£ that can be viewed as a region 
where the errors in the lattice gas model are dominate. 
Yet, at the relevant larger scales (far above the small 
scale ^), the physics of the quantum lattice gas model is 
indistinguishable from that predicted by the relativistic 
quantum field theory representation of Dirac fields. 

Yet, in the context of the toy model, we should use ([M)) 
to calculate the vacuum energy associated with a spin-i 
fcrmion. When we integrate over all space to determine 
the total density contained in a Dirac field we have 



Ptot 



(Pk 

c2(27r)3 



27r2c£4 



V(pc)2 + (mc2)2 



(25a) 



dike) (/c^) Vsin^ k£ + cos'^ H(25b) 
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a sub-Planckian length scale.' 



wavelength (1.616 x 10 m) 

FIG. 1: Log-log plot of (|24|l for mass (red dots) and mo- 
mentum (blue dots) in GeV of a single proton versus its 
wavelength with the small scale set to the Planck length, 
£ = £p = 1.616 X 10"^^m. The straight lines are the de 
Broglie relation of quantum mechanics, p = h/X (blue dashed 
line), and the invariant mass of special relativity, mo = 0.511 
MeV (red dashed line). Respectively, the slopes are —1 and 
for the standard theories. The intersection of the mass and 
momentum lines occurs at the Compton wavelength of the 
Dirac particle. The two insets are linear plots in the extreme 
ultraviolet region, ^/4 < A < 201, where the lattice theory 
departs from quantum mechanics and special relativity. 



where the quantity e = moc£/h is small when nio is much 
less than the mass, nio <C h{£'^ /t)~^ . If we take e = 
to model a massless relativistic particle, then performing 
the integration of (|25bp yields 




FIG. 2; Comparison of theoretical predictions of the vacuum en- 
ergy density as a function of wave number cutoff times the Planck 
length, kci. The dashed curve is the standard quantum field theory 
prediction and the solid curve is the quantum lattice gas prediction 



l|26|l . where 



.theory 



at fee 



4.08557. The theoretical predic- 



tion at kct ~ 1 is 10^'^^ times too large for both curves, whereas for 
the solid curve kct ~ 4 is very close to the experimentally observed 
value. 



Alternatively, one can choose the Planck scale to be 
smaller than the grid scale, £p <^ £ <^ \. In this case, 
we still have p^^°"='(4.08557M) = 0. There exists a real- 
valued number C < 4.08557 for which pt^^{Ckct} = , 
although we have not predicted this number and thus 
do not address the fine-tunning problem. Considera- 
tions regarding an additional fundamental length scale, 
in addition to the Planck scale, have recently appeared 
in Ref. |15| . including references therein. 



{2k£ sin k£ - {k'^£^ - 2) cos k£) 



(26) 



which can be either positive, zero, or negative depending 
sensitively on the value of the wave number cutoff kc as 
well as on the value of the small scale £. 

If one takes the quantum lattice g simplis- 
tic representation of Planckian scale physics (choosing 
t = t-p) where ([M)) is interpreted as physical behavior 
instead of numerical grid error, then its prediction of an 
allowable small value of the vacuum energy of a mass- 
less spinor field may be viewed as a new physical mech- 
anism. That is, if we use ([M)) to calculate the vacuum 
energy associated with a spin-i fermion, then the toy 
model can avoid the cosmological constant problem — it 
is not necessarily 10^^^ times too large. Taking kc as a 
parameter in (|26| . at kc = (4. 08557. ..)/£, we find that 
P™c°'^ = 0, as shown in Fig. [5] The experimentally ob- 
served value. Pill' = 9.9 x 10~2'^kg/m^ (such as obtained 
by the Wilkinson Microwave Anisotropy Probe) is ob- 
tained at a slightly smaller wave number cutoff, corre- 
sponding to a grid scale about four times smaller than £, 



VI. CONCLUSION 

We revisited the quantum lattice gas model with a uni- 
tary evolution operator UsUc applied at a small scale 
£ that advances a Dirac spinor field, represented on a 
grid, forward by a small time scale increment t. We de- 
rived the conditions for which the generator of evolution 
is the Dirac Hamiltonian, UsUc = e -"(-'^-p-+'^- ™o)^ -y^e 
quantified the error of the quantum lattice gas model as 
a departure from standard quantum mechanical behavior 
for the particle momentum going as p = {h/£) sm(2'K£/X) 
and as its departure from a relativistically invariant par- 
ticle mass going as m = mo cos(27r^/A). In this regard, 
the quantum lattice gas model ([2|) is numerically accu- 
rate only to scales > 20£, even though it retains covariant 



^ Our simplistic estimate is further simplified by not considering 
the inflationary epoch of space under an extremely high vacuum 
energy density, just the dynamics of a fermionic field when the 
spacetime is flat with a small positive cosmological constant when 
its discreteness below the Planck scale becomes relevant. 



6 



behavior down for scales > i. Yet, the numerical error 
of the model can be taken as a good feature — providing 
a mechanism for a small positive cosmological constant. 

There have been a number of theoretical attempts em- 
ploying, for example, supersymmetry Il2| , string the- 
ory [13[, and the anthropic principle 14| to bridge the 
known chasm between the large quantum field theory 
prediction of p^f^ ~ 10"'^^*'eV* for the latc-timc, zero- 
temperature vacuum energy density of "empty" space 



and the observed value of p° 



10 eV associated 



with a small positive cosmological constant. The toy 
model presented herein gives a value of fcc for which the 
vacuum energy vanishes. This is the root of the equation 
2kcismkc£—{k'^£'^ — 2) cos kci—2 = 0. However, the theo- 
retical considerations presented above do not tell us why 
the wave number cutoff should be fined tuned so that 
P^^aT'^ — Pvac Nevertheless, the quantum lattice gas 



model appears to be one potential route to reconcile a dis- 
cretized quantum field theory, at least a version modified 
at a small scale by ([24| . with the well accepted exper- 
imental observation of a positive cosmological constant 
by employing a plausible wave number cutoff parame- 
ter that corresponds to a fundamental grid scale. In a 
subsequent paper, we will numerically evaluate the novel 
unitary collision operator (jl8b[) employed in a quantum 
algorithmic simulation of the dynamical behavior of a 
system of Dirac particles. 
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